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Abstract Some fractional and anomalous diffusions are driven by 
equations involving fractional derivatives in both time and space. 
Such diffusions are processes with randomly varying times. In rep- 
resenting the solutions to those equations, the explicit laws of cer- 
tain stable processes turn out to be fundamental. This paper di- 
rects one's efforts towards the explicit representation of solutions 
to fractional and anomalous diffusions related to Sturm-Liouville 
problems of fractional order associated to fractional power func- 
tion spaces. Furthermore, we study a new version of the Bochner's 
subordination rule and we establish some connections between sub- 
ordination and space-fractional operator. 
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1 Introduction and main results 

In recent years, many researchers have shown their interest in fractional and anoma- 
lous diffusions. The term fractional is achieved by replacing standard derivatives 
w.r.t. time t with fractional derivatives, for instance, those of Riemann-Liouville or 
Dzhrbashyan-Caputo. Anomalous diffusion occurs, according to most of the signifi- 
cant literature, when the mean square displacement (or time-dependent variance) is 
stretched by some index, say a ^ 1 or, in other words proportional to a power a of 
time, for instance t a . Such anomalous feature can be found in transport phenomena 
in complex systems, e.g. in random fractal structures (see Giona and Roman [22]). 

Fractional diffusions have been studied by several authors. Wyss [60], Schneider 
and Wyss [56] and later Hilfer [26] studied the solutions to the heat-type fractional 
diffusion equation and presented such solutions in terms of Fox's functions. For the 
same equation, up to some scaling constant, Beghin and Orsingher [6]; Orsingher 
and Beghin [49] represented the solutions by means of stable densities and found 
the explicit representations only in some cases. Different boundary value problems 
have also been studied by Metzler and Klafter [45] ; Beghin and Orsingher [5] . In the 
papers by Mainardi et al. [36, 38, 39] the authors presented the solutions to space- 
time fractional equations by means of Wright functions or Mellin-Barnes integral 
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representations, that is Fox's functions. See also Mainardi et al. [40] for a review on 
the Mainardi- Wright function. 

For a general operator si acting on space, several results can also be listed. Nig- 
matullin [48] gave a physical interpretation when srf is the generator of a Markov pro- 
cess whereas Kochubei [30, 31] first introduced a mathematical approach. Zaslavsky 
[61] introduced the fractional kinetic equation for Hamiltonian chaos. Baeumer and 
Meerschaert [1] studied the problem when szf is an infinitely divisible generator on a 
finite dimensional space. For a short survey of these results see Nane [47]. 

In general, the stochastic solutions to fractional diffusion equations can be realized 
through subordination. Indeed, for a guiding process X(t) with generator stf we have 
that X(V(t)) is governed by d^u = &/u where the process V(t), t > is an inverse or 
hitting time process to a /3-stable subordinator (see Baeumer and Meerschaert [1]). 
The time-fractional derivative comes from the fact that X(V(t)) can be viewed as a 
scaling limit of continuous time random walk where the iid jumps are separated by 
iid power law waiting times (see Meerschaert and Scheffler [43] ; Metzler and Klafter 
[46]; Roman and Alemany [53]). Results on the subordination principle for fractional 
evolution equations can also be found in Bazhlekova [4]; Bochner [11]. 

Anomalous diffusions can also be carried out by considering a fractional operator 
acting on the space. The problem of finding a suitable representation for a fractional 
power of a given operator A defined in a Banach space X has a long history. The first 
definitions of fractional power of the Laplace operator were introduced by Bochner 
[11]; Feller [20]. For a closed linear operator A, the fractional operator (— A) a has 
been investigated by many researchers, see e.g. Balakrishnan [3]; Hovel and Westphal 
[27] ; Komatsu [32] ; Krasnosel'skii and Sobolevskii [33] ; Watanabe [59] . Although the 
methods presented differ, each of those papers was primarily based on the integral 
representation 

-(-A) 8 /= H \ a - x (XI-A)- l Afd\ 

n Jo 

for a well defined / and < 5R{a} < 1 under 

(i) A € p(A) (the resolvent set of A) for all A > 0; , , 

(ti) ||A(AJ - A)~ 1 \\ < M < oo for all A > 0. [ } 

Different definitions can be also given by means of hyper singular integrals, see e.g. 
Samko [54]. 

In both cases, time and space fractional equations, the explicit representations 
of the law of stable processes and, those of the corresponding inverse processes, are 
fundamental in finding explicit solutions to fractional problems. 

In this paper we study time and space fractional problems involving the operator 
Q* (see formula (3.1)) which is the adjoint of an infinitesimal generator of non-negative 
diffusions Q. In particular, for tv(x) — x 1 ^ 1 , the second order differential operator 

^=^Tt|-^" 7+1 |-. 7 = ±1, M>0 
■y z m(x) ox ox 

is the operator governing two related diffusions, the squared Bessel process G M (t), 
t > (for 7 = +1) and its inverse process E^t), t > (also known as reciprocal 
process of G^, for 7 = —1). Due to the fact that P{E^{t) < x} = P{G li {x) > t} we 
refer to as the inverse of G M . 

For such processes we study the governing equations where the time derivative is 
replaced by its fractional counterpart and find solutions in bounded and unbounded 
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domains. For the time-fractional equations on bounded domains we study Sturm- 
Liouville problems of fractional order associated with fractional power function spaces. 
A complete orthogonal set of eigenfunctions (w.r.t. the weight function to(x)) arises 
naturally as solutions of the second-order differential equations involving Q under 
appropriate boundary conditions and therefore, we obtain Sturm-Liouville boundary- 
value problems associated to the killed semigroups of G M and E^. 

The fractional power of Q* (for 7 = +1, that is the governing operator of the 
squared Bessel process G^) is also examined. We find the explicit representation of 
— (— Q*) u for v £ (0, 1) in terms of Riemann-Liouville derivatives and we discuss the 
properties of the corresponding subordinated process. Thus, we obtain a represen- 
tation of the power v of the composition of non commuting operators (formula (3.1) 
below) . 

All fractional problems investigated in this work have stochastic solutions with 
randomly varying times which are subordinators. Such subordinators are denoted by 
Z? t and ££. From the fact that P{£\ < x} = P{$j» > t} we say that 2% is the inverse 
of ff t which is a positively skewed stable process with non-negative increments and 
therefore non-decreasing paths. This means that Z v t = inf{x > : 9f x ^ (0, i)} can 
be regarded as hitting time. We find that, for v = 1/n, n G N, 

^'^^(^(...^H 1 /") ...)), t>0 (1.2) 

and 

tf= [G fil (G fl2 (...G fin (t v /u)...))r, t>0 (1.3) 

for suitable choices of fi — (pi, . . . , fj, n ). Furthermore, we show that the compositions 
(1.2) and (1.3) hold for all /1 G ^™+i(n!) where 

2?l (q) = L G R!^ : p = H v = ( Vl , . . . , v n ) G N™, f[v 3 = g\ 

I/71 I In 

with m,K,g£ N. This result permit us to explicitly write the laws of fj t and £ t 
and therefore the laws of the processes subordinated by them. In particular, we obtain 
useful representations of solutions to fractional equations involving general operators 
but representing anomalous diffusions realized through subordination. 

The main results of this work are collected in Section 3. We present auxiliary 
results and proofs in the remaining sections of the paper. 

2 Introductory remarks and notations 

We first introduce the following notation: 

• s„ is the law of the symmetric stable process S% , 

• h v is the law of the stable subordinator jjj , 

• l v is the law of Z, v t which is the inverse to Sj" , 

• <?m = 9]i i s the law of the squared Bessel process starting from zero G^(t), 

• e li = g^ 1 is the law of E^t) which is the (reciprocal of ) inverse to G M (f), 

• H™' n is the Fox's function, 
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and W" are the left and right Weyl's derivatives, 



and a/ a ' _ s „ are the left and right Riemann-Liouville derivatives, 



dx 



d{-x)° 



c Df is the Dzhrbashyan-Caputo fractional derivative. 



We now introduce fractional derivatives and recall their connection with stable 
densities. The a-stable process Sf ,e , t > 0, with law s 6 v = s^(x,t), x £ M, t > has 
characteristic function 



Sexp 



exp 



i6— tan 

lei 



~2~) 



(2.1) 



with a € (0,1) U (1,2] and 6 G [-1,1] (see Zolotarev [63]). If = 0, then we have 
a symmetric process with E exp (ij3S") = exp (— i|/3| Q ), a G (0,2]. For the sake of 
simplicity we will write S" instead of S^'° . Moreover, we will refer to ff t = S"' 1 , t > 
as the totally (positively) skewed process which is also named stable subordinator. 
For n — 1 < a < n, n E N, according to Kilbas et al. [29]; Samko et al. [55], we define 



and 



(W°f)(x) 



(-1)" d n 

r(n — a) dx n 

1 d n 

T(n — a) dx n J_ 



(s - x) 



n—a — 1 



f(s) ds, 



x € 



(x-s) 



n— a— 1 



f(s) ds, x G 



(2.2) 



(2.3) 



which are the right and left Weyl's derivatives by means of which we write the gov- 
erning equation of S"' e , t > 0, given by 



dt 



(x,t) = e D" x \S a (x,t), x G E, t > 



where 



kl 



2 cos an/2 



kW? + (1- k)W1 



(2.4) 



(2.5) 



and < k = k(9) < 1 (see e.g. Benson et al. [7]; Chaves [13]; Mainardi et al. 
[36]; Orsingher and Zhao [50]). For = (that is k = «(0) = 1/2) we obtain the 
Riesz derivative 

which is the governing operator of the symmetric process S", t > 0. The Riemann- 
Liouville derivatives 

"'"' f : (x) = (W*f)(x), xeR + (2.7) 



d(-x) a 



and 



dx a 



(x - s y 



7(s) ds, x G 



(2.8) 



1 d n 
r (n — a) dx™ Jg 

are dehned by restricting the Weyl's derivatives (2.2) and (2.3) to (0, +oo). For a G N, 
the fractional derivatives above become ordinary derivatives and 



dx 



- = (-l) c 



d a 
d{-x) a ' 



(2.9) 
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We also introduce the Dzhrbashyan-Caputo fractional derivative 



"~D?f(t) 



1 



r (n - a) Jo 



(t-s) 



n-a-1 d n f 

ds n 



(s) ds 



(2.10) 



defined for n — 1 < a < n, n 6 N, which is related to (2.8) as follows (see Gorenflo 
and Mainardi [24] and Kilbas et al. [29]) 



s w(*) = S(*)-E 



dt k 



(*) 



fc=0 



j-k—a 



t=0+ 



r(fc-a + l)' 



(2.11) 



From the relation 

Pr{£? < = Pr{S% > t], (2.12) 

according to [1; 16; 42; 44], we define the inverse process £j, £ > with law l v = 
l v (x, t), x, t > 0. As already mentioned before, Sj", t > is the ^-stable subordinator, 
v G (0, 1) with law, say, /ij, = fej,(a;, t), x,t > 0. The process , i > is a process with 
non-negative, independent and homogeneous increments (see Bertoin [8]) whereas, 
the inverse process t > has non- negative, non-stationary and non-independent 
increments (see Meerschaert and Scheffler [43] ) . Stable subordinators and their inverse 
processes are characterized by the following Laplace transforms: 



Eexp {-XSjt) = exp (-tA" 



£7exp(-A£ t ") = E u (-M") 



and 



C[h v {x, -)](A) = x^E^i-Xx"), C[l v (x, -)](A) = A^ 1 exp {-x\ v ) 



where the entire function 

E a3 {z) = Y, 



k>0 



T(ak 



-, z G C, ft{a} > 0, G 



(2.13) 
(2.14) 

(2.15) 



is the generalized Mittag-Leffler function for which 

A a-/J 



/■CO 

y o e - A v- i ^(-M a )dz v , + r 



, 5R{A} > |c| 1/Q , 5R{c} > 



and E a (z) = E a ^{z) is the Mittag-Leffler function. From (2.13) we immediately verify 
that the law h v satisfies the fractional equation — t) = ~g^ph v (x,t) whereas, 

for the law of £,", from (2.14) we have that mslv{x,t) - 
laws can not be represented in closed form. In this paper we write 



^lu(x,t). Such density 



h v {x,t) = -^Hil 



X 




t l/u 


(0,1) 



» lu(x,t) = ^Hl'i 



(1 - v, v) 
(0,1) 



(2-16) 

for x, t > 0, v G (0, 1) in terms of H-functions as can be obtained by considering 
(A. 12) and the Mellin transforms (see [16; 36; 37]) 



M[h v {-,t)]{ v ) =r 



1 — 7] 

V 



tV 1 r (n) vto- 1 ) 



Further representations of h v and are given in terms of the Wright function 

„fc 



k>0 



MT(ak + f3) 



zeC, 5R{a} > -1, P G 
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by considering that 

l v {x,t) = ^W- v , 1 - v (--) 
and (see D'Ovidio [17]) xh v (x,t) = tl v (t,x). 



3 Fractional and Anomalous diffusions 

Standard diffusion has the mean squared displacement (or time-dependent variance) 
which is linear in time. Anomalous diffusion is usually met in disordered or fractal 
media (see e.g. Giona and Roman [22]) and represents a phenomenon for which the 
mean squared displacement is no longer linear but proportional to a power a of time 
with a / 1. Thus we have superdiffusion (a > 1) or subdiffusion (a < 1) in which 
diffusion occurs faster or slower than normal diffusion (see e.g. Uchaikin [58]). 

3.1 Fractional evolution equations 

We begin our analysis by studying anomalous diffusions on Q a — (0, a), a > 0, whose 
governing equations involve the operator 

Q* = 1 9 x 7A1 " 1+1 9 1 (3.1) 
7 2 dx dx vo{x) 

(tu(x) = x 1 ^ 1 will play the role of weight function further on) for 7 = ±1 and fi > 0. 
The solutions to the ordinary problem 

dt 

subject to the initial data uq = f can be written as u(x, t) — T\(t)f(x) where T\(t) = 
exp—tg*. For the subordination principle (see for example Bazhlekova [4]; Bochner 
[11]) we can write the solutions to the fractional problem 

is defined in (2.8)) subject to the initial condition uq — f by considering the 
convolution operator 



1 f 00 
T v {t) = ~ J o dsH 



1,0 
1.1 



(0,1) 



where the Fox's function H\ ,{ 1 is the law of the inverse process introduced in the 
previous section. In order to explicitly write the convolution T u (t)f(x), for 7 ^ 0, 
/i > 0, we introduce the functions 

ff2(»,t)=sign(7)i0j(|) and = g^x,^) (3.2) 



where 



W = ^YUN e ~ z ' z>0 '7>o, /i>o 



is the well known generalized gamma density or Wcibull distribution if 7 € N. 







Theorem 1. For 7 7^ 0, /i > 0, v 6 (0,1], the solution to the fractional p.d.e. 

^§^ = g*u2'", x&n 00 ,t>o 

m2^(x,0) = 5(x) 



(3.3) 



is given &?/ 

u™(x,t)= I gl(x,s)l v (s,t)ds (3.4) 



where g 1 and l v are defined in (3.2) and (2.16) respectively. 



From the fact that 



lim £? = t 



where t is the elementary subordinator (see Bertoin [8]) we can write 

lim l v {s, t) = S(s — t) 

and thus uj'^ = g 1 are the solutions to (3.3) for v — 1 and 7 7^ 0, [i > 0. In this case 
the time-fractional derivative becomes the ordinary derivative d/dt. For 7 = ±1 and 
/i > 0, we obtain that 

~gl{x,t)=E*8{G lx (t)), 

~g- 1 (x,t) = E*5(E^t)), 

ul<»(x,t) = E*6(G fl (L» t )), 

u~^(x,t)^E x S(E^m)) 

where E x 5(X t ) = 6 * fx t (x). The process G M is a non-negative diffusion satisfying 
the stochastic differential equation 

dG^t) =fxdt+ 2y f G f Jt) dBi{t) (3.5) 

where Bi(t), t > is a Brownian motion with variance t/2. The reciprocal gamma law 
e p = ff^ 1 ' M > 0, represents the 1-dimensional marginal law of the process satisfying 
the stochastic equation 



dE^t) = - (e^) - dt +] j 2| ^ (t i )|2 dB 2 (t) (3.6) 

where B^{t), t > is a standard Brownian motion (see e.g. Bibby et al. [9], Peskir 
[51]). Due to the global Lipschitz condition on both coefficients, the stochastic equa- 
tion (3.6) has a unique solution which is a strong Markov process. The process E^ 
also appears by considering the integral of a geometric Brownian motion with drift 
/x, that is 

1 r°o 

(3.7) 



1 f°° 

~E M '= / exp (2B(s) - 2/is) ds 

2 Jo 



see Dufresne [19]; Pollack and Siegmund [52]. For 7 = 2, the operator Q* becomes 

1 /2 

the governing operator of a 2 /^-dimensional Bessel process i?2^ = . 

For the processes G^ and E^ introduced above there exist a couple of interesting 
properties. In particular, we have that 

Pr{G^x)>t}=Pr{E^(t)<x} 
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and therefore we refer to as the inverse to G^, whereas from the fact that 

E,(t) la = w l/G^t) 
the process will be also termed reciprocal process of G^. 

Remark 1. We observe that ff t has non-negative increments and therefore, from 
(2.12), the inverse process £^ can be regarded as an Hitting time. This does not hold 
for E^ being G M a diffusion driven by (3.5). 



The operator Q* is the adjoint of the infinitesimal generator 

-, 7 = ±l,/i>0 (3.8) 



1 9 yw-r+l 9 



j 2 w(x)dx dx' 

which is a second order differential operator driving the squared Bessel process 
if 7 = +1 and the inverse process E^ if 7 = —1. The Sturm-Liouville eigenvalue 
problem (see Bochner [10]; Courant and Hilbert [14]) associated with (3.8) leads to 
the differential equation (see Lemma 2 below) 

Gi> Ki = -(Ki/2) 2 i; Ki . (3.9) 

The eigenfunctions 

i, Ki {x) = xW-riJu-i (n lX ^ 2 ) (3.10) 

corresponding to different eigenvalues are orthogonal with respect to the weight func- 
tion vo{x) = x 1 ^ 1 in the sense that 

^ (D ^ (?) w{x)dx = °' if * ^ 3 (3 - n) 

on the domain £l a = (0, a), a > 0. The eigenvalues are written in terms of {KijieN 
which are the zeros of the Bessel function of the first kind (see [35, p. 102]) 

^ = r4a + Hi) ' w< °°" |argz|<7r - (3 - l2) 

The set of eigenfunctions {tp Ki }ie^ is complete and therefore a piecewise smooth 
function can be represented by a generalized Fourier series expansion of (3.10). In 
particular, on the finite domain Oi = (0, 1), we study the solution to 

c Dlm™ =g*ml>i*, xefli,t>0, 

mWfaO) =mo(x), m eC(ni) (3.13) 
ml^\x,t) = 0, x e dQi, t > 0, 

with v € (0, 1], 7 ^ 0, fi > and arrive at the next result. 

Theorem 2. TTie solution to the problem (3.13) can be written as follows 

m™(x, t) = n{x) £ c n E v (-(n n /2) 2 t") ff^ (3.14) 

where Xo(x) = x 1 ^ 1 is the weight function, 

mo(x)ip Kn (x) dx, n=l,2,... (3.15) 



'Hi 

and E v (z) = E Vt \{z) is the Mittag-Leffler function (2.15). 
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This result can be easily extended to the case in which the bounded domain is 
O a = (0, a) for a > and we can write 

m£'"(M) =E x mo (G^))M^<T SiA G,)) (3.16) 

and 

m-^(x,t) = E x m (E^))l {£ , <Tna(E ^ (3.17) 

where T D (X) = inf{i > : X t £ D} is an exit time and E x (j){X t ) = <f> * /x t (a:)- We 
observe that l(z,"<T n (X)) = l(t<Tn (x(L")))- F° r a — > oo we obtain the solutions to 
the problem (3.3). Indeed, we have that l(Lj-<T n=0 (G^)) = 1 an d l(ij , <T fioo (E (1 )) = 1- 
Furthermore, formula (3.11) can be rewritten for a — > oo as 

/>00 

/ V5 Ki (a:)V' Kj (x)ro(a;)da; = 8{n % - 
Jo 

which leads to the Hankel transforms 

/•OO j* oo 

(Hf)(p)= xJ„(px)f(x)dx and f(x)= pJ u {px){H f){p)dp 
Jo Jo 

of a well-defined function /. 



3.2 Fractional powers of operators 

From the Cauchy integral 

27T? 7 r (z - Z ) 

we can define an algebraic isomorphism such that a function of a linear bounded 
operator A is defined as 

f(A) = ±Jf(z)R(z,A)dz 

where R(z,A) — (zl — A)^ 1 is the resolvent operator (under conditions (1.1)). In 
general, for a closed linear operator A in a Banach space X, the definition of A a for 
a complex a could be given by means of the Dunford integral 

A a =— [ ^^-A)- 1 dC = -^^ [°° X a (X + A)- 1 dX 

where T encircles the spectrum cr(A) counterclockwise avoiding the negative real axis 
and C° takes the principal branch (see e.g. [3; 27; 32; 33; 59] and the references 
therein). For such operators the expected property A a A 13 = A a+ @ holds true. A 
well-known example is the fractional Laplace operator which can be also defined (in 
the space of Fourier transforms) as 

A Q / 2 u (x) = -^/ e -**||e|r.F[u]«)d£, xel" (3.18) 

The stochastic solution to the Cauchy problem involving the fractional Laplace op- 
erator (3.18) is given by the process B(f)") (B is a Brownian motion driven by the 
self-adjoint Laplace operator A and f)" is a stable subordinator) and has been first 
investigated by Bochner [11]; Feller [20]. In the one-dimensional case, the operator 
(3.18) becomes the Riesz operator (2.6) for which the representation (2.5) given by 
means of the right and left Weyl's derivatives holds as well. 
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In this section we study a fractional power of (3.1) which is the power of the 
composition of non commuting operators and obtain the explicit representation 

A f = -(-£*)" / = (x 1 - /)) (3-19) 

defined on the positive real line fioo = (0, +oo) for v £ (0,1). This representation 
involves the Riemann-Liouville fractional derivatives (2.7) and (2.8) which replace the 
Weyl's derivatives (2.2) and (2.3) for functions defined on the positive real line. As we 
can check, from (2.9), the fractional operator (3.19) for v — 1 becomes the operator 
(3.1). We show that, for v £ (0, 1), the process G^Sj") is the stochastic solution to 
the Cauchy problem involving the operator (3.19) where G M is driven by Q* . First we 
state the following result. 

Theorem 3. Let us consider the process G^{t), t > 0, satisfying the stochastic equa- 
tion (3.5). For /jl > 0, W £ (0, 1], the 1— dimensional density law g^ of the process 
solves the fractional p.d.e. on U,^ — (0, +oo) 

subject to the initial condition g^ix, 0) = 8{x). 

The stochastic solution to (3.20) is the process G^ which does not depend on the 
fractional index v £ (0, 1]. Furthermore, we notice that the space operator appearing 
in (3.20) is different from A. 

We proceed our analysis by considering the process F"'^ = Srf fi , t > with law 

/•OO 

f Vi p(x ) t)= h u (x,s)l p {s,t)ds, x> 0, t > 0, i/,/3 £ (0,1) (3.21) 
Jo 

which has been thoroughly studied by several authors, see e.g. [16; 18; 28; 36; 38]. If 
v = j3, then the law (3.21) takes the form f VlJJ (x, t) = i -1 f„(t~ 1 x) where 

Mx)= 1 ^T 18 ! 11 !! *v > *>M>0, «/e(0,l) (3.22) 

7T 1 + 2X V COS TTU + X iv 

and F"' u l = t x i^/ 2 ^t> * > °> (which means that F"' v £ P x ) where 3 = 1,2 
are independent stable subordinators and the ratio iSj" / 2^)" 1S independent oit (see 
[12; 16; 34; 62]). The density law /„ arises in many important contexts, we refer to 
the paper by James [28] and the references therein for details. 

Lemma 1. The governing equation of the process F" , t > 0, with density law 
(3.21), is written as 

w + ^)^ Ax ' t) = 6(x) w^~y x ^°> t>0 ( 3 - 23 ) 

with fv^tdClcojt) — and f„ i; g(iE, 0) = 5(x) or, by considering (2.11), 

{ cd * + £0 fi/Ax ' t] = °' x - o> 1 > (3 - 24) 

with fu.p{x,Q) = S(x). 
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Proof. From the Laplace transforms (2.13) and (2.14) we have that 

A) = J e- xt (E e-« F °'") dt = A"- 1 ^ + £ v )-\ 
Let us consider the equation (3.24). From the fact that 

C [ c D?f] (A) = A^£[/](A) - A^V(0 + ), p G (0, 1) 
(which comes from (2.11)) we obtain 

A^*(e,A)-A^ 1 +r*(C,A) =0 
which concludes the proof. 

We state the main result of this section concerning the operator (3.19). 
Theorem 4. For x G fi^ = (0, +oo), t > 0, fi > and f3, v G (0, 1] we have that 
i) the density law 0„(a;, t) = J °° gj^ix, s) h v {s, t) ds solves the fractional p.d.e. 



□ 



dt 



(x,t) = AQ v Jx,t), 



(3.25) 



ii) the density law g"'^(x,t) — J °° g\{x, s) fu,p{s,t) ds solves the fractional p. d. 



dtP 



(x,t) = Ag v /(x,t). 



Furthermore, 



1 



v v ^(t t) - 1 r v 'P ( x 



(3.26) 



(3.27) 



where 



(i,£); (1=1); (m,o) 

( M ,l); (1,|); (1,1) 



a: > 



is a Fox's function defined in (A. 13). 
From Theorem 4 we have that 



»(x,t) = E*5(G^ t )), 
S f(x,t) = E*S(G,(F^)) 



and thus Gy,{F" ), t > is the stochastic solution to (3.26) whereas G^{F^' L ) = 
G^Sjf), t > represents the stochastic solution to (3.25). This is because of the fact 
that Sj^i =' ft", t > 0, being £* =' t the elementary subordinator and l v (s,t) —> 
5(t— s) for v — > 1. Furthermore, from the formulas (A. 12) and (A. 11), by taking into 
account (4.30), we immediately have that 



p X 



CM) 

CM); 

(i,£); GM) 
(m,i); 



, x > 0, v G (0,1), 
, a: > 0, i/ G (0, 1) 
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and 



rl,0 



CM) 



5G > 0. 



For /3 = v, the equation (3.26) takes the form 

1 M—hK^w^^"^^) <3 ' 28) 



<9f 

which differs form (3.20) and represents the governing equation of the process 

Gfj, (t x i^/ai^i), « >0 

where jSj", t > 0, j = 1, 2 are independent stable subordinators. Here, the stochastic 
solution to (3.28) depends on v only by means of the ratio which possesses 

distribution (3.22). For (3 — 1, equation (3.26) becomes (3.25) and writes 



at 



(fl^ = gj^' 1 ) with stochastic solution given by G^(fi"), t > 0. Finally, for ^ = 1 
in (3.26), we reobtain the governing equation of the process G^L"), t > already 
investigated in Theorem 1. 

Remark 2. We observe that, for /i£N, 

G M (itf) = ||B(ioHI| 2 =E Wrf 

3=1 

where B(t) = (i_B(t), . . . , n B(t)), t > and jB(t), j — 1, . . . , n are independent Brow- 

nian motions. From the Bochner's subordination rule we get 3 -.B(j5£) = jS 2 " which 
are symmetric stable processes with E expi^jSf" — exp— i|£| 2l/ for all j = l,,..,n. 
Thus, by considering n independent stable processes jSf" and S 2ly = (iSf v , . . . , nS" 2 "), 



i > 0, we obtain that G^Sft) = ||S 



Remark 3. For B(t) € R 3 , the three dimensional Bessel process represents a ra- 
dial diffusion on a homogeneous ball. The subordinated squared Bessel process can 
be therefore regarded as a radial diffusion on a non-homogeneous ball, with fractal 
structure for instance. This interpretation is due to the fact that the random time Sj" 
has non-negative increments and therefore, non-decreasing paths. Furthermore, the 
subordinated process G^fy") speed up as 9f t increases. For v —> 1, Sf t — > t a.s. and 
GpiSy^ ) becomes standard diffusion because of the linear growing of time. 

Remark 4. From (4.31) we obtain that 



E 



G^W) cxt^, r >0, v,0e (0,1]. 



Moreover, — > £f for ^ — > 1 and therefore we have that G„(£j ) is a subdiffusion 
whereas, from the fact that Fj — > ijj for /3 — > 1 we get the superdiffusion G^fjJ'). 
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3.3 Explicit representations of solutions via Mellin convolu- 
tions 

The laws of the processes S) v t and £^ can be written in terms of H-functions as pointed 
out in Section 2. Alternative expressions can be given in terms of the Wright function 
but only for v = 1/2, 1/3 we obtain a closed form of the density laws and therefore 
of the solutions investigated in the previous section. 

We give an explicit representation of the solutions presented so far by exploiting 
the Mellin convolutions of generalized gamma functions. The simplest convolutions 
we deal with are written below: for x, t > and j ^ 0, fii, fj,% > 0, 

9l *ffu 7 M) = , ' ' \ 7 — ir" (3.30) 

where £?(•, •) is the Beta function (see e.g. Gradshteyn and Ryzhik [25, formula 8.384]) 
and, 

2 N (^ 7 A 7 ) 

where K a is the modified Besscl function of the second kind. In particular 



v(\ n I- a (z) - I a (z) 

Ka{z) — — ; , a not integer (3.32) 

2 sin an 

(see [25, formula 8.485]) where 

( z /2) a+2k 

W = E fcl nL + fc + iv N<TO ' l ar s 2 l <7r ( 3 - 33 ) 

k>a ' v ' 
is the modified Bessel function of the first kind (see [25, formula 8.445]). 
Definition 1. For — oo < a < b < oo, we define the space 

M b a = {/ : R + h. C | x^fix) e L 1 (R+),Wr ] € M b a } 
where U b a = {C : CeC,o< 9?{C} < 6}. 

Definition 2. Let us consider the function Vy : C 4 K. We define the class of 
one- dimensional processes 

P Q = |y(t), < > : 3§Ci' s.t. Ely^)/*!"" 1 = TVfa). Vt? G §} , aeE. 

Remark 5. We notice that S% eP„ ^ S% l = t x l v S\. 

Definition 3. We define the class of functions 

F« = {f\Y ~ f,Y GP Q } 

where Y ~ / means that the process Y possesses the density law f . 

Remark 6. We remark that F a C Mj. 

We point out that for a composition involving the processes Y a . = \X aj ] a where 
X a . are independent processes such that X a . € P Q for all j = 1, 2, . . . n, we have that 
Yctj G Pi, Vj which implies that 

r CT1 (y CT2 (. ..Y an (t) . . .)) ^ Y ai (t 1 ^)Ya 2 (t 1/n ) ■ ■ ■ ^„(i 1/n ) (3.34) 
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for all possible permutations of {<Xj}, j = 1,2, ... ,n. This can be easily carried out 
by observing that Vy a . iv) — ^Px a . (v a — a + 1). Indeed, 

and therefore, for i =/= j, 

E[Y aj (Y ai (t))]"- 1 = V x „. (m-a+ l)E[Y ai (t)]"" 1 = E[Y a] (t 1 / 2 ) Y a% (t 1 ' 2 )]^ 1 . 

From the fact that Pi 9 ~ gl G Fi and 

P 7 9 [G^h- ~qI gj G F 7 

Pi 9 [G^)] 1 ^ ~ 5 2 ^ fl^Fx 
we can define the following Mellin convolution of gj G Mf^_ „. 

Definition 4. Jbr 7^0, /x = . . . , fx n ) E R™ we define the Mellin convolution 

gl>* n (x,t)=gl 1 *...*g1Jx,t) (3.35) 

with Mellin transform (see (A.l)) 

M[g^ n (-,t)](r 1 )=f[M{gl.(-,t 1 / n )]( V ) = f" 1 f[ r(( ^ 1)/ ^ + ^ ) (3.36) 
3=1 3=1 1 

where n E and a = 1 — minj{7/Xj}. 

We notice that f x « *f Y e = /yfl */x» is the law of X a ■ Y 13 if X G V a and Y G Fp 
whereas the well-known Fourier convolution fx * fy is the law of X + Y. Also, we 
introduce the sets 



S?k) = { p E R n + : p = -, v = (m, . . .,v n ) E N n , J2vj = ? \ (3.37) 



and 



'2( Q )={peRl : p=-,«=K,.. A )6N", 11^' = ^ (3 - 38) 

K 3=1 



with m,K,g E N. For 7 = 1,2 and a fixed = (/ii, . . . , G =5^™(<j), we have that 

9h*~-* 9l n (x, t)=g 1 g ^*---* g1 an (x, t) 

for all — (8 ai , . . . ,8 an ) E =5^™(<r) and all permutations of {crj}, j = 1, 2, ,..,n. 
This fact follows easily from the semigroup property (*-commutativity) of the law g 1 
which will be shown in ii), Lemma 5 below. We observe that H = |^™| < |N| is 
the cardinality of &™ , thus &™ is a finite set. Furthermore, Vg E N and a fixed 
fj, E Z?™(q), we have that 

M[g^(;t)](v) = M[ 9 r n (;tm), V0G^( e ) (3.39) 
whereas, for \x E S fi ^ l {<;) and 7 = 1,2, we have that 

•%r n M)KO = W(->*)](0. e (3.40) 

where we used, the familiar notation, /*" = f fll * ■ ■ ■ * f^ n . The symbols M. and T 
stand for the Mellin and Fourier transforms. 
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Theorem 5. Let us consider v = l/(n + 1), n 6 N and fi = . . . , /x n ). 

ij For the stable subordinator ft" , t > 0, the following equivalence in law holds 
true 

E fll (E^(...E fln (( V t)^) ■■■)), t>0, M G^" +1 (n!) (3.41) 
where the process Eu(t), t > 0, satisfies the SDE (3.6). 
ii) For the inverse process 2% , t > 0, the following equivalence in law holds true 

[G, 1 (G, 2 {...G, n {t u /v)---))W t>0, /xe^ +1 (n!). (3.42) 
where the process G^t), t > 0, satisfies the SDE (3.5). 

From Theorem 5 and formulas (3.30), (3.31) we can explicitly write g^'* 9 . Further- 
more, this representation holds in the general set ^™ +1 (n!). Indeed, for v = l/(n+l), 
n £ N, the stochastic solution to (3.25) is given by 

G fi (E IXl (E IM2 (...E lJtn ((vt) 1 /") ...)), t>0, ( ftl ... )Mn )e^ +1 (n!). 

Moreover, the representation (3.42) turns out to be useful in representing the solution 
to the problems (3.3) and (3.13). A natural extension follows for the problem (3.26). 

From E^cb) = \Efjt, t > 0, c> and formula (3.41) we obtain that 

S% ^ v-^E^E^. ..E^{t 1/y ) ■ ■ ■)), t > 
and thus, for a± — l/(ni + 1), a-i = l/(ri2 + 1) and 

Mi = (Aii,i,...,Mi, ni ) e ^(m!), m 2 = (M2,i,---,M2,nJ e ^" 2 2 +1 (n 2 !) 

we have that 

F t aua2 '= v-^E^E^i- ■ ■ E„ 1>ni (G Ml 2 (G M2 2 (. . . G M2 ,„ 2 (t7„) • • •))) ■••))■ 

(3.43) 

The fact that E^ e Pi and G M <E Pi means that 

E^G^it)) 1 ^ G^E^t)) 

or equivalently 

From this we can write the law (3.21) in terms of the convolutions (3.30) and (3.31). 
Corollary 1. For v = l/(n+ 1), n £ N, t/ie stochastic solution to (3.26) is given by 

G^), t>0 

where G M has law and F v '^ has density which can be represented by means of the 
Mellin convolutions (3.30) and (3.31) as formula (3.43) entails. 

Remark 7. For 7 ^ 0, /1 > and v — 1/5, the stochastic solution to (3.3) can be 
written as follows 



G M (G Ml (G P2 (G AI3 (G A14 (t 1 / 5 /5))))) 



1/7 



t > 
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(G M , G Ml , . . . , G M4 independent squared Bessel processes) where 

p = (jh, to, Ha, IH) € ^ 5 4 ( 4! )- 
For /x = (3, 2, 2, 2)/5, we have that 

1/5 V > ; f2 / 5 J o J Q s 2/5 5^ s y l) ^53/2 f 1/2 

where 

2r(|)r(f) 



C 



5 3/ 2 [*r(§)rr(§) 

For further configurations of fi 7 see Remark 9. 

Remark 8. The relation between stable densities and higher order equations has 
been investigated by many authors (see for example Baeumer et al. [2]; DeBlassie 
[15]; D'Ovidio [17]; D'Ovidio and Orsingher [18]). In [17] we have shown that the law 

li of solves the higher-order equation 

1 ' dx n dt' 

4 Auxiliary Results and Proofs 

4.1 The operators Q and Q* 

The operators we deal with are given by 



and 



Vfl.v 1 slsw A|1 ' /3£D(r) (41) 



x 1 -! f d 2 d \ 



^- 7+ i . , flED{ g) (4.2) 



7 2 to(x) dx \ dx 

where to(x) — x 7At_1 . We shall refer to Q* as the adjoint of Q. Indeed, as a straight- 
forward check shows, we have that Q* to fi = to Q fi and the Lagrange's identity 

f-iGh- hG* fi = Q (4.3) 

immediately follows. Thus, by observing that 

D(G*) = {/ e Mi : / = w/i, A e £>(£?)}, 

(see Definition 1) we obtain that 

(Gfi, h) = (fx, G*f 2 ), V/i G £>(£) and V/ 2 G 
Lemma 2. 27ie following hold true: 
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i) For the operator appearing in (4.2) we have that 

^ k = (k/2) 2 Vk (4.4) 

where 

iP K (x) = xi^-^K^ (o l/2 j , k>0,x> 0,7=^0 (4.5) 
and K a is the Macdonald's function (3.32). 
ii) For the operator appearing in (4.2) we have that 

Q4> K = -{K/2fi> K (4.6) 

where 

$ K (x) =x^ x -^J il - l (k^/ 2 ) , k>0 iX > 0,7^0 (4.7) 
and J a is the Bessel junction of the first kind (3.12). 

Proof. We first recall some properties of the Macdonald's function (3.32): we will use 



d a 
—K a {z) = -K a -x{z) - -K a (z). 
dz z 



(4.8) 



(see [35, p. 110]); the functions K a and I a are two linearly independent solutions of 
the Bessel equation 



^d 2 Z QL {^X^ dZ a (^X^j 



dx 2 



dx 



x Z a (x) = 



(4.9) 



whereas, the functions J a and Y a (see [35] for definition) are linearly independent 
solutions to 



2 d 2 Z a (x) dZ a {x) 2 

+ x Z a (x) = v 



(4.10) 



dx 2 dx 

(see [35, pp. 105 - 110]). 

By performing the first and the second derivative with respect to x of the function 

ipn = ipK,(%) we obtain 



^ =1(1 + ^(1-/^7/2 



2x 



-K_ 



1- A* 

KX"*/ 2 



¥ - - IH (2 ~^ -¥- - -¥x x ^ )K -» 



and 



K 

KX^I 2 



2x 



By keeping in mind the operator Q , from the fact that 

2^2 



(4.11) 
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the relation (4.4) is obtained. The equation (4.11) can be rewritten as 

x 2 + ( 7 /x - 7 + 1) x - 7 2 («/2) 2 ^ V« = (4.12) 

which is related to the formula (4.9) whereas, a slightly modified version of (4.12), 
which is 

x 2 i>'L + {l^-l + l)x4,' K + 1 2 (k/2) 2 4> K = 0, (4.13) 

is related to the formula (4.10). The equation (4.0) can be written as formula (4.13) 
and therefore, after some algebra, from (4.10), we have at once that 

as announced in the statement of the Lemma. □ 
Formula (4.13) can be put into the Sturm-Liouville form as follows 

(^- 7+ V«)' + 7 2 W2) 2 to(a:) 4 = 0. (4.14) 

According to the Sturm-Liouville theory ([14]) and formula (4.14), we obtain an or- 
thogonal system {V'ftilisN such that 

e^« i = -(«i/2) 2 ^ (tl , (4.15) 

where Ki are the zeros of J a and Q is a Hermitian linear operator whose eigenfunctions 
are orthogonal w.r.t. the weight function Vo(x) — x 7 ^" 1 . Indeed, from the fact that 

/ x J v (nix) J v (kjx) dx = 0, if i^j 
Jo 

(see [35]) we get that 



ii 



i> Ki (x)ip Kj (x)to(x)dx = 0, if i ^ j. (4.16) 



4.2 Proof of Theorem 1 : time-fractional diffusions in one- 
dimensional half-space. 

For v = 1 the density law (3.4) becomes the law of G M , uj'^ = whose Mellin 
transform is written as 

/ \ — 

y t (v)=M{gl(;t)]( v )=r^l—+^ j/eHJV (4.17) 

We perform the time derivative of (4.17) and obtain 

d _ , , 77 — 1 _, / 77 — 1 \ '7-T-i 



9i 7 V 7 

7/-1 / 77 - 7 - 1 + 7/A ^ /'»? - 7 - 1 



7 V 7 / V 7 

=-2-(»y - - 7 - 1 + i^t(v - 7) 

7 
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= \(V - 1)(V - 7)*t('7 - 7) + -^(V ~ 1)(7M - l)*t(»? - 7) 



1 



7 7W 



7 Z 



— x 2 ' 1 — a 1 
dx dx M 



fo) - 



7 



dx X 9p 



(ry). 



From the fact that £ Mi and according to the properties (A. 8), (A. 2) and (A. 3), 
the inverse Mellin transform yields the claimed result. We give an alternative proof 
by exploiting the Laplace transform technique. The Laplace transform of gl(x, t), 
x, t > 0, can be evaluated by recalling that (see [25, formula 3.478]) 



x v 1 exp {— f3x p — 72; p \ dx = - ( ~ 

p Vp 



(2^) 



(4.18) 



where p, 7, /?, z/ > and if^ is the modified Bessel function. Thus, we obtain 



c[m{x,-)){\) =2 



x 2 



(H-i)-i 



r(/i)A 



(2A 1 /V/ 2 )=2^gl / (A)^(x;2A 1 /2 ) 



where /(A) = A^ and ip R (x) = ^>(:r; k) is that in Lemma 2. By considering that 
g2(x,t) = to(x)k^(x,t) and ^* ro(a;)fc^ = to(x) Qkl we get that 

L[G* gy(x, -)](A) =2^M /(A) Q ^(x; 2X 1 / 2 ) = A -)](A) 

where in the last formula we used the result (4.4). From the fact that 



C 



(\) = \C[gZ(x,-)}(\), x>0 



we obtain the claimed result for v = 1. 



Now, we consider v £ (0, 1). From the Laplace transform 
C[l v (x, 0](A) = A"" 1 exp(-xX v ) 
(see formula (2.14)) we obtain that 



/WOr, 01(A) 



ro(x) A 



r(M) A 
to (a;) 



r(/i) 



where = Vfo «) is th at in (4.7) with k = 2X V I 2 and /(A) = A 1, (" +1 )/ 2-1 . Thus, 

in the right-hand side of (3.3) we obtain 

£[<?* u™(x, 0](A) = f ^ /(A) ro(x) ^(x; 2X»' 2 ) = 2^|/(A) £ 1>(x; 2X»' 2 ) 

where we have used the fact that Q* to / — to £/ /. Finally, from (4.4), we obtain 

C[G* u™(x, 0](A) = X" C[ul>"(x, OKA). (4.19) 
We note that |u2' M (-,t)| — Be~ qot for some B, go > as a function of t and thus, 



(A)=A^[^(x,0](A), 



(4.20) 



see [29, Lemma 2.14]. By comparing (4.19) with (4.20) the result follows. 
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4.3 Proof of Theorem 2 : regular Sturm-Liouville problems. 

From the fact that Q* m^'^x, i) = \o(x)Q m2' M (x,i), the problem (3.13) reduces to 

J^m™ = Sm™, m^(dn 1 ,t) = 0, m™(x,0)=m (x)/tt>(x) (4.21) 

where m1' ,1 {x,t) = to(x) ra2' M (x,i) and Vo(x) — x 7/i_1 . Furthermore, from Lemma 2, 
we have that Q ip Ki — — (Ki/2) 2 ip Ki where 

i> Ki {x) = x^ 1 -^.]^ (W /2 ) (4.22) 

and Ki, i € N are the zeros of J a . Formula (4.16) leads to the orthonormal system 
{^ Ki (a;)/||V' Ki ||^; ieN} where ||/||^, = (/, f) m is the norm associated to the inner 
product (4.16) with respect to the weight function ro(x). Thus, 

oo 

L 2 (R) =0^„ 
n=i 

where W n is the space of eigenfunctions associated with the eigenvalue X n = (k„/2) 2 
and we obtain that 

m^{x,t) = f^c n (t,X n )fsM- (4.23) 
where ||V> K J| ro = J'n-i{ K n) I s/ 7 ! ( see J e -g- [35, p. 130]). From (4.21) we have that 



which holds term by term. From the fact that 

G^n n {x) = -X n ip Kn (x) 
where A n = (k„/2) 2 (see (4.6)), formula (4.24) lead to the fractional equation 

c Dt c n (t, X n j = — X n c n (t, A n ) 

and thus, we obtain 

c n {t,X n ) = c n -E v {-X n t v ) (4.25) 

(the Mittag-Lefher is an eigenfunction of the Dzhrbashyan-Caputo fractional deriva- 
tive) where c„ must be determined by taking into account the initial data. In partic- 
ular, 

c n = (wo/w, </> K „)ro = / m (x)'ip Kn (x) dx. 
Formula (4.23) solves (4.21) and we obtain 



fhZ'^x, t) = {m™(x, •), U; t)) = J2 c n E v (-( K „/2) 2 f ) (4.26) 
We have to observe that m^'^x, i) = tt>(x)m 7 ,'' i (x, t) for the proof to be completed. 
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4.4 Proof of Theorem 3 

For fi > 0, a G (0, 1) 

3ScHg°s.t. (TI^ a kl)(v)=0, (4.27) 

(see the Appendix A) where gl(x, t) =Xo(x) k^{x,t) and kj(x,t) = \j\/T(fi) exp(— (x/<) 7 )/i 7AI . 
Indeed, being 

( I o- a f)( x ) = Y^jJ x (s-xr- 1 f(s)ds, x>0 

we obtain 

(4r a *i(-> *))(*) ^-^^x,*) 

and (4.27) immediately follows. We restrict ourselves to the case v G (0, 1). From the 
formula (A. 9), we obtain 



- v ) Jo 



r "" r^-i^-^^-^-ppj^*))) & 



rfo) r^-n-i & 



OO 



r(r?) r^+M-i) 
"r(r?-^ r^ + M-i-^io 

r(T,- V )r(r, + n-l-v)J 9,{x,t)dx 
T(rf) r(?7 + ^- 1) ! 



r(r) — v) T(rj + fj, — 1 — v)~ 
The x-Mellin transform of both members of (3.20) writes 

where M[gl(; t)](r}) = t^rfa + /z- l)/r(/z), 77 G Hf%. Thus, we have that 

opMMMv) = r{r) _ v) m * 

because of the fact that j^t ' 1 = T(0)/T(j3 - a)^- 1 (see e.g. [55, Property 2.5]) 
which concludes the proof. 

4.5 Proof of Theorem 4 

We proceed as follows: first of all we find out the Mellin transform of the fractional 
operator acting on space 

MM = -^{ X ^ 1+V d^^" fM ))' VE[0A) (428) 
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for a well-defined function / G M^°, a <E K (and for which (4.27) holds, that is 
(T/o± a /)(??) = 0) and second of all we prove ii) by exploiting the Mellin technique 
and then i) as a particular case of ii). We also consider / € D(Q*). Let us write 

= W, \ and *„(//) - 



r(i-r?) " wy r(»i + /* — 1 — i/) 

From (A. 10) we have that 

x*i- 1 Af(x,t)dx = -$ u (r l ) J x *-»-i x »- 1 +»——^(xi-i*f( x ,t))dx 

roc av 

= -$ u (v) J x^- l 2^(x l -»f(x,t))dx. 

From (A. 9) we obtain 

x n - 1 Af{x,t)dx = -$ u {rj)^ v {rj) / x r > + > i - u - 1 x 1 -' i f(x, t)dx 

Jo 

= -® v {r,)* v {ri)M[f{;t)]{r 1 -v). 
Thus, by collecting all pieces together we have that 

«*) - £lTr\ w,M,-*Y (4») 

Now, we consider the x-Mellin transform 

X fe/ 0, *)] (v) =M [gl(; 1)] ( V ) x M [f^ (•, t)] (v) 



where the fact that h u € F„ and lp € Fjy^ leads to 

M[f^(',i)](l)=^[i(',l)](l)^[W^)] 
and, from the formulas (2.17) we obtain 



v 



M \k u 'H- t)] M - T ^ + t i -V r (V) r (V + 1 ) t ^P new 1 (4 30) 

where a = max{0, 1 — /i}, > 0. Now, we show that 

^M[g^(-:*)]W = M[A?/(-,t)](v) (4-31) 

by taking into account the formula (4.29). The right-hand side of the formula (4.31) 
can be written as 



r^(. au^ - _ r(i -?? + !/) r( ?? + M -i) 

T(l- V ) r(ri + n-v--\ 

r(^ + /i-i)(V)r(V) r(a=i) 



r( M ) ^r(i-r?) r (2^/3-^ + 1) 
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r^ + M-i) r(^) r(^ + i) 



r( M ) v r(i-ri)r(*=±p-i3 + i) 



=Mlx v > f} (- 1)1 fn) T ( V v 1 P +1 ) 



From the fact that 



d^n=lp = r(^ + i) ^ g _ 
dtp r(*=±p-p + i) 

(see [55, Property 2.5]) formula (4.31) immediately follows and this prove ii). 



For j3 = 1, the formula (4.30) takes the form 



M[ U;t)](v) = 



ur(jA)r(i-ri) 



r 



where a = max{0, 1 — /i}. This is (for (i = 1) because of the fact that fj^i a =' Sj", 

t > 0, being a =' i the elementary subordinator, see e.g. Bertoin [8]. Thus, form 
(4.29), the Mellin transform of both members of (3.25) becomes 



where 



^•■')io-)- ^ ( 7';:V 



1 - 77 + 1/ 



?7 — 1 \ T(?7 + jU - 1 - v) 

v ) ttvT[h) F(l — T] + v) 
T(rj + fi — 1 — v) 

ttv r(/i) r(i — rj + v) 



r 



1 — r/^ 9 n=l 
St* " 



*l[G^(-)]fa)=M$ 



By collecting all pieces together we obtain the result claimed in i). 
From (4.30) and by direct inspection of (A. 12) we arrive at 

V 

where G v /[x) = ^/(x, 1). Thus, 

1 



(i-M); (i-f>f); (m,o) 
(i (o,i) 



2.1 

3,3 



(m-i,i); I); (o,i) 

(1,1); (1,§); ( M ,0) 
( M ,l); (1,±); (1,1) 



where we used, for c = 1, the property of the H functions 



Tj~m,n 
p,q 



(b 3 ,l3j);j = l,..,q 



1 



Trm,n 
P,1 



(a % + ca t , aj)j=i,.. lP 
{bj+<^j,Pj)j=\,.., q 



(4.32) 



for all c e R (see Mathai and Saxena [41]). From (A. 2), by observing that 
1 



M 



( 



( V )=M{G»/(-)}(r ] )t 



we obtain the claimed result. 
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4.6 Proof of Theorem 5 

Hereafter, we extend the result given in [16] (Lemmas 3 and 4) and show how the 
Mellin convolution turns out to be useful in order to explicitly write the distributions 
of both stable subordinators and their inverse processes. Let us consider the time- 
stretching functions ip m (s) = ms 1 ^ 771 , s £ (0, oo), m £ N and </? m such that ip m — (p^ 1 
(the inverse function of <p m ). 

Lemma 3. [16, Lemma 2] The Mellin convolution e™(x,tp n +i(t)) where fij — jv, 

for j = 1,2, ... ,n is the density law of a v-stable subordinator {Ljj: , t > 0} with 
v = l/(n + 1). neN. Thus, we have 

h v (x,t) = e*J l (x,ip n+1 (t)), x,t>0. 



We recall that = g^ 1 is the 1-dimensional law of E^. 

Lemma 4. [16, Lemma 3] The Mellin convolution g^ +1 ^* n (x,ip n+ i{t)) where fij = 
j v, j = 1, 2, . . . ,n and v = l/(n+ 1), n £ N, is the density law of a v-inverse process 
{i[ , t > 0}. Thus, we have 

l„(x,t) = gfr + V'* n (x,ip n+1 (t)), x,t>0. 

We observe that is the 1-dimensional law of Gf^ {t n+1 ) = G l [ 1 {t 1 / V ). 

The following facts will be useful later on. 
Lemma 5. For = gj(x,t), x £ M.+, t > 0, /i > the following hold: 

i) -k-commutativity: g^\ * gf 2 = gj 2 2 * gj\ for all 71,72 ^ 0. 

ii) *-commutativity: g 1 ^ * g l2 2 — g"® * for all 71,72 7^ 0. Furthermore, for 
7 = 1,2, 

* 9fj, 2 — 9^+^12- 

Hi) (*,*)-distributivity: when -k- and*- commutativity hold, we have 

9^ * {g^ 2 *9l 3 ) = (ah *9j 2 ) * (six *£# 3 )- 

iv) for /Ui,/i 2 ,c £ N 

9l., 2 = *£ =1 & = ™ d 9l = ^. M2 . c±1 (4-33) 

where #~ =1 / = /1 * f 2 *, ■ ■ ■ , *fn- 

Proof. The point i) comes directly from the formula (3.34) and the fact that gfy £ Fj, 
V7j ^ 0, and fij > 0, j = 1,2. We show that ii) holds. For 7 = 1, V* > 0, gj 
is the gamma density with Laplace transform C[g\{-, i)](A) = 1/(1 + Xt)^ and the 
statement follows easily. This is a well-known result. The case 7 = 2 is considered 
in Shiga and Watanabe [57] being g 2 , the semigroup for a Bessel process R 2fl = G\^ 
where G M satisfies the stochastic equation (3.5). The result in Hi) can be obtained 
by considering, Vt > 0, the independent r.v.'s G'JfAt), j = 1,2,3 with densities 
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gfi. = gfcixtt)) j = 1,2,3, x £ R + . From the fact that gfy G F 1; j = 1,2,3 we have 



law 



that Gl\{Gf 2 {t)) = GJl(Gl\(t)), that is, Vj, are commutative under *. For this 
reason and the *-commutativity, Vf > 0, we can write 

= G l 2 ( G U t )) + G lMM- 

In the last calculation we have used the fact that 



E 



cxp(-A[G Ml ( S ) + G M2 ( S )]) 



s = X t 



= E 



exp(-AG Ml+M2 (s)) 



s = X t 



The same result can be achieved for 7 = 2. In order to prove iv) we proceed as 



follows: first of all we observe that ii) implies g 1 ^. 
of all we consider that 



^2 



lflii- Second 



whereas 



*j"c=i *ji=i 9^2 ^ *j c =i9mj±3. - fl^.^a 
and this concludes the proof. 
Proposition 1. The following holds true 

gr n {^) = 9log^-l\x^l V Mj e /x, i = 1,2,. 

w/iere 



□ 



(4.34) 



/1 MM) 

/or : [0, +00) M-[0, +00), i = 1,2. 



fi(x,s) f 2 (s,t) ds 



Proof. Fix n = 3. Vt > 0, it is enough to consider the r.v.'s GJ, G 7 and their 
density laws gj, g^ where g 7 (x,i) = g^(x,t 1 ^' 1 ) or equivalently G^(f) G^(t 7 ). 
In this setting, we have that X{t) = G'^ 1 (GJ l2 (GJ l3 (t))) can be written as X(t) == 
G^(Gi 2 (Gi 3 (P))) thank to the fact that (Gj) 1 l = Gj r Thus, we can write 

9T 3 (x,t) = & ° ^,)(^* 7 ) = ^ (si, *<4)(m 7 ). 

Thanks to the ★-commutativity we have that g^ 2 * = g 1 ^ ★ and also that 

flr 3 ^'*) = ~9l f&a)^* 7 ) = & ° (fii *<4)(M 7 ). 
By considering n processes, the formula (4.34) immediately appears. □ 
Remark 9. For v = 1/5, we have that 

ii/sfo *) = 5 * 1/5 )' Mi = (1/5, 2/5, 3/5, 4/5) (4.35) 

and, from the Proposition 1 and the Lemma 5, 

5,*4 _ ~5 1,*3 

9(1/5,2/5,3/5,4/5) — #1/5 ° %/5,3/5,2/5) 



=01/6 ° 



(9 1/5 * 9l/ 5 * 9l/5 * ffl/s) * (fl/5 * 9l/S * 5l/5) * 32/5 
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=9l/5 ° 
=01/5 ° 



/ l l l l \ / !>* 2 M 2 i,* 2 \ 

l9l/5 * 9l/5 * 01/5 * 01/5J * 19(2/5,1/5) * 9(2/5,1/5) * 0(2/5,1/5)-' 

( 1,*3 \ 
1 ^(2/5,1/5,1/5); 



k= 

_~5 1,*3 / r5i\ _ 5,*4 

— 01/5 ° 024/5, 1/5, 1/5V X > <7 — 9(24/5,1/5,1/5,1/5) 

where 

M 2 = (24/5,1/5,1/5,1/5) 
and /x l5 /x 2 G (24). Finally, we obtain 

5-*4 / ,\ 5,*2 5,*2 / ,\ 

0(24/5,1/5,1/5,1/5)^' l ) —9(24/5,1/5) * 9(1/5,1/5) \ X i 

From (3.31) the corresponding integral representation emerges. 

Theorem 6. Fix /x G ^™ (g). TTiera #b*™ = g*'*" /or a// 19 G ^ (e). 

Froo/. Fix K,g£t We have = g^-k . . . * fl£ n , M = (Ml, • • • Mn) e From 
(3.38) we can write /1 = ...,//„). Let us first consider n = 2. We recall that 

9(/i* 2 /i 2 ) = 9(jLi2 /ii) f rom ^ ne *-commutativity. Thus, form the properties i) and m) of 
the Lemma 5, we have that 

1,*2 _ Ai 1,*2 _ Ai 1,*2 

9(i Al ,i A2 ) ~ *ii=i 9 ( i,i A2 ) - *j 1= i9 ( i A2! i) 

_ ^Ai ¥ A2 l,* 2 _ l,* 2 
- *j 1= i *, 2= i9 ( a^A) - */ 1= i *j 2= i 0a (m) - 

For n G N, i?o = ^(1) . . . , 1) G , we can write 

l,*n _ l,*n _ Ai $ A™ 

9^ -%£!,...,#„) - *j'i=i ••• *i»=i 5 A> • 

We shall refer to #0 as the O-configuration. We first observe that 

n 

l,*n _ * An _ _l,*n „ _ TT n ■ 

j=l 

or equivalently 

l,*n _ Ai vJin _ _ 

9/x - *j 1= i • ■ • *i„=i 9^ - *ii=i%e/£ t ,i,....i) - 9i(g i ,/j i ,i,...,i) 

ft = p/Aii f° r au * = 1, 2, . . . , n. The last identity comes from the ★-commutativity. 
By exploiting the *-commutativity and the *-commutativity we have that gji* n = 
gg'* n for all 6 such that 



1 

t™9 0=-( & ji i ,i), ft=e/n 



K 

Sj=l 



where dim(l) = n - |i| - 1, fa = (£ 8l) . ■ -,fa w ) € N^ 1 , Sj G i, j = 1,2,... |i| and 
|i| < 71 is the cardinality of i. A further configuration is given by 6 = (gi,fa)/K where 
|i| = n — 1. In this case, 6 G t3?™(g) is obtainable by n! permutation of the elements 
of fx. In a more general setting, for a = {ot\ ■ /3, a 2 , ... , a n ) G N™, /3 G N, c G N the 
following rules hold 

9a = *i=i9 ( ; i)Q2 ,... )Q „) ( 4 -36) 



2G 



(see iv), Lemma 5) and, for c > 1, 



^/3±c 1,-A-n 



^/3±c 1,-kn 1,-kn 
' j=l"(ot < , 1 ,a<r 2 ,...,a<r n ) ~ 9(a ai - 



(4.37) 



for all permutation of {<Xj}, J = 1, 2, 



We recall that 



for a,P,c G N. By making use of the properties i), ii) and Hi) of the Lemma 5 we 
can obtain all possible configurations of i? G starting from the O-configuration 
i?o- All different configurations of $ are included in ^"(g). From (4.37), for all 
■d = . . . , # n ) G ^"(q) we have that = 6- This concludes the proof. □ 

By collecting all pieces together we obtain the claimed result. 



A Fox functions and Mellin transform 

The Mellin transform of / G (see Definition 1) is defined as 

/>oo 

M[f(-)]( V )= / x^ x f{x)dx, V €M b a . 



(A.l) 



Let us point out some useful operational rules that will be useful throughout the 
paper: for some -co < a < b < oo and b > 0, /, /i,/2 G M^: 



x^f{bx)dx =b~*>M [/(•)] (r?), 



* G) * 



X[x b /(-)] fa) [/(•)] (v + b), 
ds 



2 IS 



fa) =X [/i(-)] (.V) x M[f 2 (-)] fa), 



where 



I(x) 



f(s)ds, x > 0, 



(A.2) 
(A.3) 
(A.4) 
(A.5) 

(A.6) 



see e.g. Glaeske et al. [23]. The formula (A.4) is the well-known Mellin convolution 
formula which turns out to be useful in the study of the product of random variables. 

We say that / G M n if / G and is a rapidly decreasing function such that 

d k f 

3o€lsi. lim x^^-^-ix) = 0, fc = 0, l,...,n, n G N, x G K+ 

x->+oo ax 



and 



3b G K s.t. lim x b 



dx k 



(x) = 0, k = 0, 1, . . . ,n, n G N, x G 



For / G M n _i and n G N we have that 



A4 



-(■) 



(A.7) 
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T(l + ra — 77) 
T(l- V ) 



M [/(•)] (v-n). 



(A.8) 



As a generalized version of the integer derivatives (A. 7) and (A.8) we introduce the 
Mcllin transform of fractional derivatives (2.7) and (2.8) (see Kilbas et al. [29]; Samko 
et al. [55] for details). For a given / e M* and < a < 1, if ^{77} > 0, 



M 



d a f 



(■) 



_d(-x) a 
whereas, if 3?{?7} < a + 1, 



fa) 



T(r/ - a) 



M[f{-m-a) + {TllZ a f){ri) 



M 



dx a 



(■) 



(V) 



r(l + a-Tj) 



M [/(•)] (T?-a) + (T/^ a /)(»?) 



where 



and 



T(l-ry) 

1 



(A.9) 



(A.10) 



r(a) 
1 



r(a) Jo 

are the right and left fractional integrals 



(s~x) a - 1 f{s)ds, x>0, 
(x- s) a - 1 f(s)ds, x>0 



The Fox functions, also referred to as Fox's H-functions, H-functions, generalized 
Mellin-Barnes functions, or generalized Meijer's G-functions, were introduced by Fox 
[21] in 1996. Here, the Fox's H-functions will be recalled as the class of functions 
uniquely identified by their Mcllin transforms. A function / G can be written in 
terms of H-functions by observing that 



n Trm,n 



((2j , 0^)2=1,., ,p 



dx 
x 



with r\ G i[j where 
M 



4 771,71 

P,Q 



nr=i r (^+^)nr=ir(i 



rioii) 



UUm+i r(i - bj - vfr) Utn+i r(«, + mi) ' 



(A.11) 



(A.12) 



Thus, according to a standard notation, the Fox H-function can be defined as follows 

M2£{ri)x-"dq (A. 13) 



Trm.n 

p,q 



1 

2tt7 



where V(M^) is a suitable path in the complex plane C depending on the fundamental 
strip (H„) such that the integral (A. 11) converges. For an extensive discussion on this 
functions see Fox [21]; Kilbas et al. [29]; Mathai and Saxena [41]. 
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